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Abstract 



Massless diquark jet and gluon jet productions at large transverse momen- 
tum in the coherent diffractive processes at hadron colliders are studied in the 
two-gluon exchange parametrization of the Pomeron model. We use the he- 
licity amplitude method to calculate the cross section formulas. We find that 
for the light quark jet and gluon jet production the diffractive process is re- 
lated to the differential off-diagonal gluon distribution function in the proton. 
We estimate the production rate for these processes at the Fermilab Tevatron 
by approximating the off-diagonal gluon distribution function by the usual 
diagonal gluon distribution in the proton. We also use the helicity amplitude 
method to calculate the diffractive charm jet production at hadron colliders, 
by which we reproduce the leading logarithmic approximation result of this 
process we previously calculated. 

PACS number(s): 12.40.Nn, 13.85.Ni, 14.40. Gx 
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I. INTRODUCTION 



In recent years, there has been a renaissance of interest in diffractive scattering. These 
diffractive processes are described by the Regge theory in terms of the Pomeron (IP) exchange 
The Pomeron carries quantum numbers of the vacuum, so it is a colorless entity in QCD 
language, which may lead to the "rapidity gap" events in experiments. However, the nature 
of Pomeron and its reaction with hadrons remain a mystery. For a long time it had been 
understood that the dynamics of the "soft pomeron" is deeply tied to confinement. However, 
it has been realized now that how much can be learned about QCD from the wide variety 
of small- x and hard diffractive processes, which are now under study experimentally. In 
Refs. the diffractive J/ip and T production cross section have been formulated in 

photoproduction processes and in DIS processes in perturbative QCD. In the framework 
of perturbative QCD the Pomeron is represented by a pair of gluon in the color-singlet 
sate. This two-gluon exchange model can successfully describe the experimental results 
from HERA §. 

On the other hand, as we know that there exist nonfactorization effects in the hard 
diffractive processes at hadron colliders First, there is the so-called spectator effect 

||, which can change the probability of the diffractive hadron emerging from collisions intact. 
Practically, a suppression factor (or survive factor) u Sp n is used to describe this effect [ID 



Obviously, this suppression factor can not be calculated in perturbative QCD, which is now 
viewed as a nonperturbative parameter. Typically, the suppression factor Sp is determined 
to be about 0.1 at the energy scale of the Fermilab Tevatron |J. Another nonfactorization 
effect discussed in literature is associated with the coherent diffractive processes at hadron 
colliders 0, in which the whole Pomeron is induced in the hard scattering. It is proved in 
that the existence of the leading twist coherent diffractive processes is associated with a 
breakdown of the QCD factorization theorem. 

Based on the success of the two-gluon exchange parametrization of the Pomeron model 
in the description of the diffractive photoproduction processes at ep colliders J2],|3|,|5j , we may 
extend the applications of this model to calculate the diffractive processes at hadron colliders 
in perturbative QCD. Under this context, the Pomeron represented by a color-singlet two- 
gluon system emits from one hadron and interacts with another hadron in hard process, in 
which the two gluons are both involved (see Fig. 1). Therefore, these processes calculated 
in the two-gluon exchange model are just belong to the coherent diffractive processes in 
hadron collisions. Another important feature of the calculations of the diffractive processes 
in this model recently demonstrated is the sensitivity to the off-diagonal parton distribution 
function in the proton [|TT|j . 

Using this two-gluon exchange model, we have calculated the diffractive J / ip production 
12fl ) charm jet production massive muon pair and W boson productions [|14]] in hadron 



collisions. These calculations show that we can explore the off-diagonal gluon distribution 
function in low x region and study the coherent diffractive processes at hadron colliders 
through these processes. In this paper, we will calculate the light quark (massless) jet and 
gluon jet productions at large transverse momentum in the coherent diffractive processes at 
hadron colliders by using the two-gluon exchange model. As sketched in Fig. 1, the two-gluon 
system (in color- singlet) emitted from one hadron interacts with another hadron to produce 
the final state light quark (or gluon) jets. There are three different partonic processes 
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contributing to the diffractive dijet productions, gp — > qqp, qp — > qgp, and gp — > ggp as 
shown in Fig. 2, Fig.3 and Fig. 4 respectively. 

In the calculations of Refs. |I2|-|14[|, there always is a large mass scale associated with the 



production process. That is for J/tp production, m c for the charm jet production, M 2 
for the massive muon production (M 2 is the invariant mass of the muon pair) and for 
W boson production. However, in the processes we will calculate in this paper, there is no 
large mass scale. So, for the light quark jet and gluon jet production, the large transverse 
momentum is needed to guarantee the application of the perturbative QCD. Furthermore, 
the experience on the calculations of the diffractive di-quark jet photoproduction |I5| shows 
that the light quark jet production in the two-gluon exchange model has a distinctive feature 
that there is no contribution from the small l\ region < kj) in the integration of the 
amplitude over l\. So, the expansion (in terms of l\jM\) method used in Refs. fT2}-|14"|j can 
not be applied for the calculations here. In the following calculations, we will employ the 
helicity amplitude method to calculate the amplitude of the diffractive light quark jet and 
gluon jet production in hadron collisions. We will show that the production cross sections 
for these processes are related to the differential (off-diagonal) gluon distribution function in 



the proton as that in the diffractive di-quark jet photoprodution process [|T^]. (In contrast, 
we note that the cross sections of the processes calculated in Refs. ]r2]-Il are related to the 
integrated gluon distribution function in the proton). 

The diffractive production of heavy quark jet at hadron colliders has also been studied by 
using the two-gluon exchange model in Ref. |fL6[| . However, their approach is very different 



from ours []. In their calculations, they separated their diagrams into two parts, and called 
one part the coherent diffractive contribution to the heavy quark production. In our ap- 
proach, in line with the definition of Ref. , we call the process in which the whole Pomeron 
participates in the hard scattering process as the coherent diffractive process. Under this 
definition, all of the diagrams plotted in Fig. 2, Fig.3 and Fig.4 for the partonic processes 
gp — > qqp, qp — ► qgp and gp — > ggp contribute to the coherent diffractive production. 

In addition, we must emphasize that in our calculations we only discuss the contribution 
to the diffractive dijet production from the so-called "Pomeron fragmentation" region, i.e., 
M\ ~ 4/cfi, where M\ is the invariant mass of the diffractive final states and kx is the 
transverse momenta of the outgoing jet. However, in the region of M\ 3> 4k^, the diffractive 
production may be dominated by the central (for M\ system) or initial hadron fragmentation 
contributions, and in this region a resolved Pomeron model (with parton distributions for 
the Pomeron) may be more appropriate. 

The rest of the paper is organized as follows. In Sec. II, we will give the cross section 
formulas for the partonic processes of Figs. 2-4 in the leading order of perturbative QCD, 
where we employ the helicity amplitude method to calculate the amplitude for these pro- 
cesses. We will also use this method to recalculate the diffractive charm jet production 
process gp — > ccp, by which we can reproduce the leading logarithmic approximation result 
given in Ref. |i~3| . In Sec. Ill, the numerical results for the diffractive light quark jet and 
gluon jet production at the Fermilab Tevatron will be given. And the conclusions will be 
given in Sec. IV. 



1 For detailed discussions and comments, please see [13] 
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II. THE CROSS SECTION FORMULAS FOR THE PARTONIC PROCESSES 



In the leading order of perturbative QCD, the Feynman diagrams for the three different 
partonic processes are plotted in Figs. 2-4 respectively. The two-gluon system coupled to 
the proton (antiproton) in these diagrams is in a color-singlet state, which characterizes the 
diffractive processes in perturbative QCD. Due to the positive signature of these diagrams 
(color-singlet exchange), we know that the real part of the amplitude cancels out in the 
leading logarithmic approximation. To get the imaginary part of the amplitude, we must 
calculate the discontinuity represented by the crosses in each diagram of these figures. 

In our calculations, we express the formulas in terms of the Sudakov variables. That is, 
every four-momenta hi are decomposed as, 

hi = atq + fop + k iT , (1) 

where q and p are the momenta of the incident parton (quark or gluon) and the diffractive 
proton, q 2 = 0, p 2 = 0, and 2p ■ q = W 2 = s. ai and $ are the momentum fractions of q 
and p respectively, is the transverse momentum, which satisfies 

h iT ■ q = 0, k iT ■ p = 0. (2) 

All of the Sudakov variables for every momentum are determined by using the on-shell 
conditions of the momenta represented by the external lines and the crossed lines in the 
diagram. The calculations of these Sudakov variables are similar to those in the diffractive 



charm jet production process gp — > ccp [pf] . 
For the momentum u, we have 



Ml 

ol u = 0, p u = xip = — ^— , u T = t = 0, (3) 

where M\ is the invariant mass squared of the diffractive final state. For the high energy 
diffractive process, we know that M\ <C s, so we have (3 U {xjp) as a small parameter. For 
the momentum k, 

k 2 

a k (l + a k ) = --rp-, Pk = -^kPn, (4) 

where k^ is the transverse momentum of the out going quark jet (or gluon jet). 

For the loop mentum /, the Sudakov parameters are not the same for every diagrams and 
every partonic processes, so we will discuss them in the following subsections respectively. 

By the form of these Sudakov variables, the cross section formula for the partonic process 
can be formulated as 

A is the amplitude of the partonic process, and we know that the real part of the amplitude is 
zero in the leading order. So, we only need to calculate the imaginary parts of the amplitudes 
for the three partonic processes. 
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In the calculations of the amplitudes for the three partonic processes, we employ the 



Helicity Amplitude method fi7Hl9f. As a cross check, we will recalculate the diffractive 
charm jet production [131 by using this method. 

Another important check to our calculations is performed by examining the behavior 
of the amplitude under the limit of l T — ► in the integration over the loop momentum 
I. From the experience of the calculations for the diffractive charm jet production in Ref. 
|13| |, we know that the contribution to the amplitude from the individual diagram of the 
parton process has linear singularity at l T — > 0. The linear singularity is not proper in QCD 
perturbative calculations. So, the linear singularities coming from different diagrams must 
be canceled out by each other for the same order of perturbative calculations. That is to 
say we must guarantee the total sum of the contributions from all of the diagrams free of 
the linear singularity in the calculations for every partonic processes. 



A. gp — > qqp process 

The first four diagrams of Fig. 2 for this subprocess are the same as those calculated 
in the diffractive photoproduction processes But, due to the existence of gluon-gluon 

interaction vertex in QCD, in the partonic process gp — > qqp, there are additional five 
diagrams (Diags. (5)-(9)). These five diagrams are needed for complete calculations in this 
order of QCD. 

First, we give the Sudakov parameters for the loop mentum I for this process (Fig. 2), 

Oil = - — , 

2(k T J T )-l T 
A = , for Diag.l, 3, 5, 

a k s 

2(k T ,l T ) + l 2 T 



(1 + a k )s 



for Diag.2, 4, 6, 

for Diag.7, 8, 9, (6) 

o 

where (k F , It) is the 2-dimensional product of the transverse vectors k T and l F . 

The imaginary part of the amplitude A(gp — > qqp) for each diagram of Fig. 2 has the 
following general form, 

ImA = C F (7%) I x Ui(k + q)T^v 3 (u - k), (7) 

where C F is the color factor for each diagram, and is the same as that for the gp —>■ cap 
process | I3J . a is the color index of the incident gluon. V ^ represents some 7 matrices 
including one propagator. F in the integral is defined as 

F = ^-/j(x',x";l T ), (8) 

where 
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dG{x',x^]l\ 



f(x',x";l T )- ' ' , (9) 



where the function G(x', x"\ h?) is the so-called off-diagonal gluon distribution function |TT 
Here, x' and x" are the momentum fractions of the proton carried by the two gluons. It is 
expected that at small x, there is no big difference between the off-diagonal and the usual 



diagonal gluon densities pOfl . So, in the following calculations, we estimate the production 
rate by approximating the off-diagonal gluon density by the usual diagonal gluon density, 
G(x',x";Q 2 ) « xg(x,Q 2 ), where x = x%> = M|/s. 

In Ref. [13[], we calculate the amplitude (|?|) for the heavy quark diffractive production 



processes by expanding T M in terms of l\. However, in the light quark jet production 
process gp — > qqp, the expansion method is not yet valid. According to the result of the 
production cross section is proportional to the heavy quark mass. If we apply this formula 
to the light quark jet production, the cross section will be zero. That is to say, the expansion 
of the amplitude in terms of Z^, in which the large logarithmic contribution comes from the 
region of l\ <C M\, is not further suitable for the calculation of the cross section for massless 
light quark jet production. Furthermore, the experience of the calculation of the diffractive 
light quark jet photoproduction process •yp — > qqp (for Q 2 = 0) |L5| indicates that there is no 
contribution from the region of l\ < k\ in the integration of the amplitude over l\. In the 
hadroproduction process gp — > qqp, the situation is the same. So, the expansion method, 
in which l\ is taken as a small parameter, is not valid for the calculations of the massless 
quark production processes. 

In the following, we employ the helicity amplitude method [|I^] to calculate the amplitude 
Eq. ([?[). For the massless quark spinors, we define 

w±(p) = -^(l±75)«(p). (10) 
For the polarization vector of the incident gluon, which is transversely polarized, we choose, 

e± = -^(0,l,±i,0). (11) 
The helicity amplitudes for the processes in which the polarized Dirac particles are involved 



have the following general forms |19 



u±{pf)Qv T {p l ) = — - =L =—, (12) 

4^(n-pi)(n-pf) 

where n is an arbitrary massless 4-vector, which is set to be n = p in the following calcula- 
tions. Using this formula (|P2"D, the calculations of the helicity amplitude A(\(g), X(q), A(g)) 
for the diffractive process gp — > qqp is straightforward. Here A are the corresponding he- 
licities of the external gluon, quark and antiquark. In our calculations, we only take the 
leading order contribution, and neglect the higher order contributions which are proportional 

M 2 

to f3 u = because in the high energy diffractive processes we have f3 u <C 1. 

For the first four diagrams, to sum up together, the imaginary part of the amplitude 
^(±,+,-) is 
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Im^(±, +, -) = a|(l + a t )Af x / ^/(x', x»; *)(- 5n? ^ >■ 

(13) 

where | and — ^ are the color factors for Diag.1,4 and Diag.2,3 respectively, and M is defined 



as 



A' p ly (14) 

-afcfl + a k ) 



The other helicity amplitudes for the first four diagrams have the similar form as 
The amplitude expression Eq. (|13|) is the same as that for the photoproduction process 
7P — > qqp previously calculated in Refs. [pj[T5| except the difference on the color factors. 
In the diffractive photoproduction process, the color factors of these four diagrams are the 
same (they are all |), while in hadropro duct ion process the color factors are no longer 
the same for these four diagrams. It is instructive to see what is the consequence of this 
difference. We know that the amplitude of the diffractive process in Eq. (^) must be zero 
in the limit l\ — ► 0. Otherwise, this will lead to a linear singularity when we perform the 
integration of the amplitude over l\ due to existence of the factor \j{l\) 2 in the integral of 
Eq. (0) fll3| . This linear singularity is not proper in QCD calculations. So, we must first 
exam the amplitude behavior under the limit of l\ — > for all the diffractive processes in 
the calculations using the two-gluon exchange model. From Eq. (0), we can see that the 
amplitude for the diffractive photoproduction of di-quark jet process is exact zero at l\ -> 0. 
However, for the diffractive hadropro duct ion process gp — > qqp the amplitude for the first 
four diagrams is not exact zero in the limit l\ — > due to the inequality of the color factors 
between them. So, for gp — > qqp process there must be other diagrams in this order of 
perturbative QCD calculation to cancel out the linear singularity which rises from the first 
four diagrams. The last five diagrams of Fig. 2 are just for this purpose. 
For example, the contributions from Diags.5 and 8 are 

• *h / \ 9 , v . , f cL 2 It „, , it 7 o \ / 1 + cxk e ( - ± - ) • {kr — (1 + oli,)It) n 

W 58 (±, +, -) = a 2 k (l + a k )M x / j^f(x', x"; l 2 T )( — * 1 1 , 2 r J >> 

J \ l t) 4 (kr - (1 + oik) It) 

(15) 

and the contributions from Diags.6 and 9 are 

W-(±,+,-) =a 2 (l + a k W, [ e ^ ( ^;"f ) ), (16) 

J \ l t) 4 (kr - oif.iT) 



and the contribution from Diag.7 is 



r rl 2 l 1 k 

lmA 7 (±,+,-)=aUl + a k )Afx J ^/(^^^(-—^l). (17) 

From the above results, we can see that the contributions from Diags.5-9 just cancel out the 
linear singularity which rises from the first four diagrams. Their total sum from the nine 
diagrams of Fig. 2 is free of linear singularity now. 
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Finally, by adding up all of the nine diagrams of Fig. 2, the imaginary parts of the 
amplitudes for the following helicity sets are, 



where 



InU(±, +, -) = al(l + a k )N x T {±) , 

lmA(±, -, +) = a k (l + a k ) 2 M x r^, (18) 



(±) _ f d% 1 2 l ^).{k T -l T ) 

l + a k e^ ±] • (fc r - (1 + a k )h) | Qfc e <±) ; (k T - a k l T ) j ^ ig ^ 



4 (fc T - (1 + a k )l T y 4 (A; T - afc / T )2 

And the amplitudes for the other helicity sets are zero in the light quark jet production 
process. From the above results, we can see that in the integration of the amplitude the linear 
singularity from different diagrams are canceled out by each other, which will guarantee there 
is no linear singularity in the total sum. 

Another feature of the above results for the amplitudes is the relation to the differential 
off-diagonal gluon distribution function fix' , x"\ I 2 -). To get the cross section for the diffrac- 
tive process, we must perform the integration of Eq. (plj). However, as mentioned above that 
there is no big difference between the off-diagonal gluon distribution function and the usual 
gluon distribution at small x, so we can simplify the integration of fll9|) by approximating 
the differential off-diagonal gluon distribution function f{x',x";l T ) by the usual diagonal 
differential gluon distribution function f g (x] 

After integrating over the azimuth angle of It, the integration T^> will then be 



T<±> = r^I, (20) 



where 

r dl 2 T 7 2 M J_/i _ k\-l\ l + Ok , 1 k T -{l + a k )l T 
J (l T y JgK ,iT) h§ K 2 2\k T -l 2 T \ }+ 4 l 2 2\k T -{l + a k )l T \ } 

4 l 2 2\k\-a k l\\ )V 1 ' 

Comparing the above results with those of the photoproduction process •yp qqp PUTS], 
we find that the amplitude formula for the diffractive light quark jet hadroproduction pro- 
cess gp —>■ qqp is much more complicated. However, the basic structure of the amplitude, 
especially the expression for the integration I is similar to that for the photoproduction 
process. In the integration of ([H]), if l\ < k T the first term of the integration over l T will 
be zero; if l T < k T /(l + a k ) 2 the second term will be zero; if l\ < k T /a\ the third term will 
be zero. So, the dominant regions contributing to the three integration terms are l T ~ k T , 
l T ~ k T /(l + ak) 2 , and l\ ~ k T /a k respectively. Approximately, by ignoring some evolu- 
tion effects of the differential gluon distribution function f g (x;l T ) in the above dominant 
integration regions, we get the following results for the integration X, 
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1 - kl) + orb* " S )] - (22) 

Obtained the formula for the integration X, the amplitude squared for the partonic 
process gp — > qqp will be reduced to, after averaging over the spin and color degrees of 
freedom, 

Pi 3 = ^Mw|('-|)' (*>) 

And the cross section for the partonic process gp — > qqp is 

da(qp —> qqp) , f ,,, 9 „ 9 9a^7r 2 1 . 2/c£.. l9 . . 



The integral bound > above shows that the dominant contribution of the integration 
over M\ comes from the region of M\ ~ Ak\. Using Eq. (Q), this indicates that in this 
dominant region is of order of 1. So, in the integration X the differential gluon distribution 
function f g (x; Q 2 ) of the three terms can approximately take their values at the same scale 
of Q 2 = kj.. That is, the integration X is then simplified to 



10MI - 27k% 



X « — x = ' T Ux; kl). (25) 
36M? m ' TJ 



Numerical calculations show that there is little difference (within 10% for hr > 5 GeV) 
between the cross sections by using these two different parametrizations of X, Eq. Q2"2"p and 
Eq. (p5|) . So, in Sec. Ill, we use Eqs. ( ^4|) and (pop to estimate the diffractive light quark jet 
production rate at the Fermilab Tevatron. 



B. Recalculate the heavy quark jet production using the helicity amplitude method 

For a crossing check, in this subsection we will recalculate the diffractive heavy quark 
jet production at hadron colliders by using the helicity amplitude method. In Ref. jnj, we 



have calculated this process in the leading logarithmic approximation of QCD, where we 
expanded the amplitude in terms of Z 2 -. Now, if we use the helicity amplitude method, we 
donot need to use the expansion method for the T M factor in Eq. ([F|) as in ||13fl . We can 
firstly calculate the amplitude explicitly by using the helicity amplitude method. 

However, for the massive fermion, the amplitude formula is more complicated. Following 



Ref. [181, we first define the basic spinors u±(k ) as, 



u + (k ) =f6iU-(ko), 
u-(k )u-(k ) = ~(1 - 75) U (26) 

where the momenta ko and k\ satisfy the following relations, 

k -k = 0, k 1 -k 1 = -1, k ■ h = 0. (27) 
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Using Eqs. fl26|) and (p7|), we can easily find that the spinor u + (k ) satisfies 

u + (k )u+(k ) = -(I+75) fa- (28) 



Provided the basic spinors, we then express any spinors u(pi) in terms of the basic ones, (18 

u ±i P.) = ^ + '^ U f''\ (29) 

It is easily checked that these spinors satisfy Dirac's equations. Now, for the massive 
fermions, the helicity amplitudes for the processes involving Dirac particles have the fol- 
lowing general forms, 



u+(Pf)Qv-(pi) 
u-(p f )Qv + (pi) 
u + ip f )Qv + (pi) 



Tr[Q(jtj - mj)(l - 75) fa{jt s + m s )} 

^(ko-Pi)(k -Pf) 
TrlQifr - m t )(l + 75) fa{jf f + m f )} 

4\/{k -pi){k -p f ) 
Tr[Q(j>i - rrij) (1 - 75) ftt fajfo + m/)] 



u_( Pf )Qv_( Pt ) = Tr ^~ m f ^(1-75) + (30) 

4y(ko-Pi)(k -Pf) 

where and my are the masses for the momenta Pi and pf respectively, where p\ = m? , p"j = 
rrif. From the above equations, we can see that for the massless fermions (m, = m/ = 0) 
the formula Eq. ( p0|) will then turn back to the formula Eq. (|T2|). 

To calculate the imaginary part of the amplitude Eq. (|7]) for the partonic process gp — > 
ccp, a convenient choice for the momenta k and k\ is, 

k =p, h = e, (31) 

where the vector e is the polarization vector for the incident gluon defined in Eq. ([Tip. 
Using the formula Eq. (^) and the above choice for the momenta ko and k\, the helicity 
amplitudes for Eq. (|?p will then be, 

InU(±, +, -) = a\{\ + a k )M x 

Inu4(±, -, +) = a k (l + a k ) 2 U x T&\ (32) 

77777, 

W(±, +, +) = W(±, -, -) = a k (l + a k )M x — (33) 
where A/" is the same as in Eq. (|T4] ) , and the integrations and Z£ are defined as 
(±) _ .., a xr f l 2 gW-k 1 g<±) ■ (fer - f T ) 

l + a k ■ (hp - (1 + a k )l T ) + q fc e^j ■ (fc r - a k l T ) ^ ^ 



4 m 2 + (k T - (1 + « fc )/ T )2 4 m 2 + (fcr _ akh y> 
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x'-i f^k. fix > x >>. pwi _ 2 ) 1 _ 1 1 



4 ml + (fcr - (1 + « fc )/ T ) 2 4 m 2 + (fcr _ ak /r)2 J 

If we approximate the differential off-diagonal gluon distribution function f(x', x"; l T ) by the 
usual diagonal differential gluon distribution function f g (x;l T ), the above integrations will 
then be reduced to, after integrating over the azimuth angle of It, 



= Tre^ • k T X c , (36) 



where 



f £kf ( 1 2 M——-—— 1 k T- m l- l T 
C ~J {l T f h{X ' T)[ lSm T 3Qk 2 72 m 2 

1 + a fc fcj, - m\ - (1 + a k )H\ a k k\ - m\ - a\l\ ^ 

+ ~3 H2 J' 



m\ cS ;/k; 



where 

m | = A^, + raj!, = v/(m|, + l T ) 2 - k T l T , m\ = J (m T + (1 + a k ) 2 l T ) 2 - (1 + a k ) 2 k T l T 



2 - ■ + Q 2J2 )2 _ ( gj 



ra 3 = \/ (ra 



and 



c 7 rJL 18ra| 36m 2 4 m 2 + 4 m| J ' 1 j 

So, the amplitude squared for the partonic process gp — > ccp will then be reduced to, after 
averaging over the spin and color degrees of freedom, 

R 2 = \o? s {^) W|4[(l - ||)4I^| 2 + ml\m (40) 

From the above results, we can see that the integrals of Eqs. (p7|) and ([39]) are proportional 
to l// 2 - in the limit of l T — > 0. That is to say that there exist large logarithmic contributions 
from the integration region of 1/R% -C l\ -C m^ for the integration over l T as in Ref. [13 . 
So, we can expand the integrals of Eqs. (j3"7|) and (|3"§| ) in terms of Zf, to get the leading 
logarithmic contribution to the amplitude. In the limit of l\ — > 0, the parameters m\, m\ 



and raj scale as, 
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raf 



ral 



An- 

mf, 


m 2 — k T 


I 2 


m\ 


i J' 
m T 


An- 

m^ 


m 2 c — k T 


(l + a k 


m\ 


m\ 


Ad- 

m T 


m 2 — k\ 


a\l\ 


m\ 


m\ 



(41) 

ra? 
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Under this approximation, the integrations Eqs. (^) and (^) will then be related to the 
integrated gluon distribution function xg(x;Q 2 ), 

2m 2 

Ic ~ 36K0W (1 ° Mi ~ 21m T )x9(x < (42) 



Substituting the above results into Eq. (pEOD, we can then reproduce the leading logarithmic 
approximation result for the diffractive charm jet production process at hadron colliders 



which has been calculated in [I3"|. 



qp — > qgp process 



We first give the Sudakov parameters for the loop mentum I for this process (Fig. 3), 



J,rp 

aii = , 

s 



2{k T ,l T )-l 



A = L for Diag .i ; 2 , 6, 



2(k T J T )+l 2 T 
(1 + a k )s 
M\-l 2 T 



for Diag.5, 7, 8, 

for Diag.3, 4, 9. (44) 



And the imaginary part of the amplitude A(qp —>■ qgp) for each diagram of Fig. 3 has the 
following general form, 

ImA = C F (T® f f^F x Ui(u - k)T, Uj (q), (45) 

where Cp is the color factor for each diagram, a is the color index of the incident gluon. V ^ 
represents some 7 matrices including one propagator. F is the same as in Eq. (^|). 

In |^T|, we calculated the diffractive photon production process qp — > •yqp, in which 



the Feynman diagrams are similar to the first four diagrams of Fig. 3. In that paper, we 
calculate the cross section by directly squaring the partonic process amplitude. However, 
in the calculations here for the partonic process qp — ► qgp because there are additional five 
diagrams contribution, it is not convenient to directly square the amplitude. Following the 
above subsections, we calculate the amplitude by employing the helicity amplitude method. 
Furthermore, we will show that by using the helicity amplitude method we can reproduce 



the cross section formula for the diffractive photon production process fzT |. 

For the massless quark spinors, we use the definition of Eq. (H^)- For the polarization 
vector of the outgoing gluon (its momentum is k + q), following the method of Ref. ]T7| |, we 
find that it is convenient to choose 

= N e [(k*+ 4) 4A^ti,)+ i 4(H ±75)]. (46) 
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The normalization factor N e equals to 



N e = — (47) 
sJ2kl 



With this definition we can easily get the scalar products between the four-momenta 
and the polarization vector e as 

k 2 s 

e ■ p = 0, e-q = N e — - — , e ■ k T = -N e k T s, e ■ l T = -N e (k Tl l T )s. (48) 
1 + a k 

And then the helicity amplitudes have the following general forms ||19|| , 

n ±{Pf )Qu ±(Pi ) = m^tjmM , (49) 
^\j{n-pi){n-pf) 

where n is an arbitrary massless 4- vector, which is also set to be n = p. Using this formula 
(HI), we can calculate the helicity amplitude vA(Ai, A2, A3) for the diffractive process qp — > 
qgp. Here Ai represents the helicity of the incident quark; A2 and A3 represent the helicities 
of the outgoing gluon and quark respectively. 

For the first four diagrams, to sum up together, the imaginary part of the amplitude 
A{+,+,+) is 

W- 4(+; +) = al{1 + x I || /(a/) x „. 4)( 2 _ ^ -{l + ak ){^ h) 

J \W 9 36 (k T - (1 + a k )l T ) 2 

(50) 

where | and ^ are the color factors for Diags.1,4 and Diags.2,3 respectively, and M q is 
defined as 

k = 1 3s = 9in. (si) 



— 2a k k T 

The other helicity amplitudes for the first four diagrams have the similar forms as (|50D , 
W 1234 (-, -, -) = W 1234 (+, +, +), 

Inu4 1234 (+, -, +) = W 1234 (-, +, -) = Zll m ^34 (+) + +) _ (52) 



These amplitude expressions Eq. ( |50D can also serve as the calculations of the amplitude 
for the diffractive direct photon production process qp — > q^p except the difference on 
the color factorsf] In the direct photon process, the color factors for these four diagrams are 
the same (they are all |). 



2 In Ref. I^TJ, we did not employ the helicity amplitude method. If we use the amplitude expressions 
Eqs.(|50"|) and ( |5^ ) (correct the color factors) to calculate the photon production process qp — ► q^p, 
we can get the same result as that in |2l|| . This can be viewed as a cross check for the methods we 
used in the calculations. 
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As discussed in the above subsection of the calculation for the partonic processes gp — > 
qqp, we must first exam the amplitude ( [iq) behavor under the limit of l\ — > in the 
integration over I to avoid the linear singularities. From Eq. (|50|), we can see that the 
amplitude for the diffractive direct photon production process qp — > q^p is exact zero at 
l T — *■ 0. However, for the process qp — > qgp the amplitude for the first four diagrams is 
not exact zero in the limit l\ — > due to the inequality of the color factors between them. 
That is to say, the sum of the first four diagrams is not free of linear singularities. So, the 
contributions from the last five diagrams of Fig. 3 must cancel out the linear singularity 
which rises from the first four diagrams to guarantee the total sum of the amplitude from 
all of the diagrams free of linear singularities. 

Finally, by adding up all of the nine diagrams of Fig.2, the imaginary parts of the 
amplitudes are 



lmA(+, +, +) = lmA{-, -, -) = -jM q x T, 

InU(+, -, +) = InVt(- +, -) = x T, (53) 



where 



d 2 l T u , „ , 2 ,Al + a k ) 2 {k T ,l T )-{l + a k )l 2 T Ak T ,l T ) + l 2 T 



(It) 9 (k T - (1 + a k )l T ) 2 (k T + 1 



r 



1 2 



{k T ,lT)-l T , 2 (k T jT) - a k l^ 
—a k — V a k — =; =i J. (54) 

(k T - l T ) 2 (k T - a k l T ) 2 

From the above results, we can see that in the integration of the amplitude the linear 
singularity from different diagrams are canceled out by each other, which will guarantee 
there is no linear singularity in the total sum. 

Following the argument in the above subsection of the calculation for the partonic process 
gp — ► qqp, we can also approximate the differential off-diagonal gluon distribution function 
f{x',x"\lj) by the usual diagonal differential gluon distribution function f g (x;l 2 -^) in the 
integration T, by which we can further simplify the integration of T. After integrated over 
the azimuth angle of It, the integration T will then be 

dlj< 2 1 + a k 1 k T — (l + ack)l T A k\ — l\ 



Q2yJ9\~>"Tji g v 2 2\k 2 T - {I + a k )l 2 T \' ' v 2 2\k T - l T \ 

1 k T — a k l\ 

2 2\k\ - a k l T] 



,%, )]• (55) 



In the above integration, if l T < k T /(l + a k ) 2 the first term of the integration over l T will 
be zero; if l\ < k\ the second term will be zero; if l\ < k^/a 2 , the third term will be zero. 
So, the dominant regions contributing to the three integration terms are l\ ~ k\j{\ + a k ) 2 , 
l T ~ k T , and l\ ~ k\ja\ respectively. Approximately, by ignoring some evolution effects 
of the differential gluon distribution function f g (x;l T ) in the above dominant integration 
regions, we get the following results for the integration T, 

r - 4) + { ^f,(*; h- <4f.fr §)]. (56) 
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Obtained the formula for the integration T, the amplitude squared for the partonic 
process qp — > qgp will be reduced to, after averaging over the spin and color degrees of 
freedom, 



—2 _ c^(47r) 3 1 + a\ 
24 M 2 x {l + a k 



\A = o/ nrr *V m 2 - ( 57 ) 



And the cross section for the partonic process qp —>■ qgp is 
da(qp -> qgp) 



Jj. 1 1=0 ; „ 

CIC JM\>Ak T 



\ dM x dk^da k [5(a k — oti) + 5(ct^ — a 2 )l 

JMl>4kl 

T\\ (58) 



e*i? 1 + e*? 1 



96(M|) 2 1 + e* fc /TT^I 



T 

M 



x 



where a,\^ are the solutions of the following equations, 



«( 1 + «) + ^ = °- (59) 

The integral bound M\ > Ak\ in Q58D shows that the dominant contribution of the inte- 
gration over M\ comes from the region of M\ ~ ^k\. Using Eq. (^), this indicates that 
in this dominant region a k is of order of 1. So, in the integration T the differential gluon 
distribution function f g (x; Q 2 ) of the three terms can approximately take their values at the 
same scale of Q 2 = k\. That is, the integration T is then simplified to 

r = 9F /s(x; + ~ 7ak + 10c ^' (60) 



D. gp — > ggp process 

For the partonic process gp — > ggp, there are twelve diagrams in the leading order 
contributions as shown in Fig. 4. The first nine diagrams are due to the existence of the 
three-gluon interaction vertex, and the last three diagrams are due to the existence of the 
four-gluon interaction vertex. But it will be shown in the following calculations, the last 
three diagrams do not contribute under some choice of the polarizations of the three external 
gluons. 

As usual, we first give the Sudakov parameters for the loop momentum I, 



for Diag.l, 4, 6, 10, 
for Diag.2, 3, 5, 11, 

for Diag.7, 8, 9, 12. (61) 



on 

A 



\rj\ 



2(k T ,l T )-l 2 T 
a k s 

2(k T ,l T ) + l 2 
(1 + a k )s 
M\~l 2 T 
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And the imaginary part of the amplitude A(gp — > ggp) for each diagram of Fig. 4 has the 
following general form, 



ImA 



C F fabc J ^2^G(ei,e 2 ,e 3 ) x F, (62) 



where Cp is the color factor for each diagram, a, b, c are the color indexes for the incident 
gluon and the two outgoing gluons respectively, and f a bc are the antisymmetric SU(3) struc- 
ture constants. G(e\, 62, e 3 ) represents the interaction part including one propagator for the 
first nine diagrams, where e\, e%, e 3 are the polarization vectors for the incident gluon and 
the two outgoing gluons. F in the integral is the same as that in Eq.(^). 
The color factors Cp for the twelve diagrams are 

C F = ^, for Diag.l, 7, 

C F = --, for Diag.2, 

Cp = —-, for Diag.3, 6, 9, 

C F = ^, for Diag.4, 5, 8, 

C F = - } forDiag.10, 11, 12. (63) 

Following the calculation method used in the above subsections, we employ the helicity 
amplitude method to calculate the amplitude Eq.(^). For the polarization vector of the 
incident gluon, which is transversely polarized, we use the same definition in Eq. (|TTj) . For 



the two outgoing gluons, we choose their polarization vectors as [17 



! ±) = iV e P+ i) (^(lT75)+^(#+ 0(l±7fe)], 

^ = N e [{i- ft 4 T 75)+ i 4{i- m ± 75)]. (64) 



The normalization factor N e has the same form as in Eq. (^7[) . Under the above choice of 
the polarization vectors for the external gluons, we can easily find that they are satisfied the 
following equations, 

p ■ e x = p ■ e 2 = p ■ e 3 = 0. (65) 

With these relations, we can further find that the last three diagrams do not contribute to 
the partonic process gp — > ggp. 

For the first nine diagrams, there are two helicity amplitudes among the eight helicity 
amplitudes do not contribute in the context of the above choice of the polarizations of the 
external gluons, i.e., 

W(+,+,+) = lraA(-,--) = 0. (66) 

In the expression of the amplitude A(X(ei), A(e 2 ), A(e 3 )), A denote the helicities for the three 
gluons respectively. The other six helicity amplitudes are divided into the following three 
different sets, 
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Ihu4(+, -) ~ lmA(-, +, +), 
Inu4(+, -, +) ~ Im„4(-, +, -), 

lmA(+,+,-)~hnA(-,-,+). (67) 

For the first helicity amplitudes set, Inu4.(±, =)=, =f), to sum up all of the nine diagrams, 
we get 

InU(±, =F, =f) = 7VVe < 1 ±) • k T X g , (68) 



where M' is defined as 



And the integration X g is 

1 f d 2 l T 2 k 2 T + (k T ,l T ) k 2 T -(k T) l T ) 

7T J [It) ( /~„ _L /„ \J 



(It) 2 (k T + h) 2 (k T -l T ) 2 

, n , N2 fc T ~ (1 + afc)(fcr,^r) 2 k^ - a k {k T ,l T ) ^ 



(A; T -(l + a fc )Z T ) 2 (fc T - afc / T )2 
1 f d 2 l T ff , 2 w, ^(Wr)+| . {k T ,l T )-l 2 T 



-(1 + afc) »(^lMzii±^S + (jfeM^], (71) 
(k T - (1 + a k )l T ) 2 (k T -a k l T ) 2 

From the above equations, we can check that there is no linear singularity at the limit of 
l T — > in the integration of the amplitude over the loop momentum. The first term of 
the integration X g in Eq.flTOD comes from the contribution of Diag.3; the second term comes 
from Diag.4; the third term comes from Diag.6 and Diag.9; the last term comes from Diag.5 
and Diag.8. The contributions from Diags.l, 2 and 7 are canceled out by each other. From 
(ffO|), we can see that the linear singularities coming from the four terms are canceled out 
by each other. The final result for the amplitude is now free of linear singularity. We must 
emphasize here that only the total sum of the contributions from all of the diagrams is free 
of linear singularity. The separation of these diagrams will cause linear singularity. 

Using the same approximations used in the above subsections, i.e., approximating the 
differential off-diagonal gluon distribution function f(x',x";l T ) by the usual diagonal dif- 
ferential gluon distribution function f g (x; l T ), we can further simplify the integration of X g . 
After integrated over the azimuth angle of It, this integration will then be 



9 = J W 19 2 " 2\k 2 T -l 2 T \ ] ~^ L + ak h 2\k 2 T - (1 + a k )l 2 T \ 



1 k T — a k l 2 
2\k T — a k l T 



+« fc (o - o, ? 2 : 2I )]- (72) 



The above equation shows that the integration X g here has the similar behavior as that of the 
integration T of Eq.(|55|) in the last subsection. So, the three terms of the above integration 
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Ig are dominantly contributed from the integral regions of l\ as l\ ~ k\j (1 + a k ) 2 , l"r ~ ^r> 
and l T ~ k T /al respectively. Approximately, we may also ignore the evolution effects of the 
differential gluon distribution function f g (x;l T ) in the above dominant integration regions, 
and so the integration I g is reduced to 

2p = 4) " (1 + «*) 3 ^; n ft , 2 ) + aj/,(x; §)]. (73) 

ft T l^i -t- a k ) a fc 

For the second helicity amplitudes set, Ihl4.(±, =p, ±), the calculations are more compli- 
cated, and the contribution from Diag.3 is 

d 2 ln 



W 3 (±, T , ±) = -AT'(1 + a fc ) 2 1 #) 



(1 TJ 

a k k 2 T ef ] ■ (k T + l T ) + {k\ + (fcr, / T ))4 ±) • 



(fcr + I 



r 



\2 



(74) 



The contribution from Diag.4 is 

lmA\±,T,±) = M'a k {l + a k ) f ^f( x ',x";l 2 T ) 

a k k 2 T ^ ] ■ (k T - Tj) + - (fc r , /r))^ ■ fcr 
(fcr - T T ) 2 

The contributions from Diag.5 and Diag.8, to sum up together, are 

W 58 (±, T ,±) = -Af'a k (l + a k ) ( ^Lf( x ', x ";l 2 T ) 

a k k T ^ ] ■ (k T - a k l T ) + (k T - a k (k T , h))^ ■ k T 
(k T - a k l T ) 2 

The contributions from Diag.6 and Diag.9, to sum up together, are 



(75) 



(76) 



Inu4 69 (±, T ,±) =M'(l + a k ) 2 J ^-f( x \ x "-l 2 T ) 



m 

a k k 2 T ^ ] ■ $t - (1 + a k )l T ) + (k 2 - (1 + a k )(k T , l T ))^ ] • k T 



(k T -(l + a k )l T ) 2 



(77) 



The contributions from other three diagrams (Diag.l, Diag.2 and Diag.7) are canceled out 
by each other. From the above results Eqs.( |7"4"|]7"7| ), we can see that every term has linear 
singularity at the limit of l\ — > in the integration of the amplitude over l T , while their 
total sum is free of the linear singularity. 

Following the procedure as we do for the helicity amplitude A(±, =F, =f) in the above, 
we can approximate the off-diagonal gluon distribution function f(x',x";l T ) by the usual 
diagonal differential gluon distribution function f g (x; l T )- After integrating over the azimuth 
angle of It, to sum up all of Eqs. (|74] - |77|) , we get the helicity amplitude, 
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W(±, =f, ±) = M'ti{1 + a k ) 2 ef ] ■ k T l gj (78) 

where X g is the same as Eq.(^) and then Eg . (|73|) under the same approximation. 

For the third helicity amplitudes set, Inv4.(±, ±, =)=), the calculations are similar to the 
calculations of Iul4(±, =)=, ±). The contribution from Diag.3 is 



lmA\±,±,T) = -Af'a k (l + a k ) f —^f( x ', x ";l 2 T ) 

(1 + a k )k 2 T ef } • (k T + l T ) - {k 2 + (k T , l T ))^ ] ■ k T 



(k T + l T ] 



2 



(79) 



The contribution from Diag.4 is 

W 4 (±,±, T ) = M'a 2 J ^f( x ',x";l 2 T ) 

(1 + a k )k 2 T ^ ] ■ (for - l T ) - (k 2 - (k T , l T ))^ ] • k T 
{k T - h) 2 

The contributions from Diag.5 and Diag.8, to sum up together, are 

d 2 h 



(80) 



W 58 (±,±, T ) = -Af'a 2 k J ^f 2 f(x',x";l 2 T ) 



(1 TJ 

(1 + a^k^e^ • (for ~ a kh) - {k T - ak{k T , h))^ • k T 

(k T - a k l T ) 2 

The contributions from Diag.6 and Diag.9, to sum up together, are 

lmA m (±,±,T)=Af'a k (l + a k ) f 

J {L T ) 

(1 + a k )k 2 T ^ ] ■ (kr - (1 + a k )h) - (4 ~ (1 + a k )(k T , h))^ ] • jfr 

(for - (1 + « fc )f r ) 2 



(81) 



(82) 



And also, we find that the contributions from other three diagrams (Diag.l, Diag.2 and 
Diag.7) are canceled out by each other, and the total sum of Eqs. (|79"1-|32]) is free of the linear 
singularity. If we approximate the off-diagonal gluon distribution function f(x',x";l T ) by 
the usual diagonal differential gluon distribution function f g (x;l T ), and integrate over the 
azimuth angle of It, their sum will lead to a similar result as in Eq. (|T8"D , 



Inu4(±, ±, t) = M'Tial^ ■ k T l g . (83) 



By summing up all of the helicity amplitudes Eqs.(^), (|78|) and |83]), we will get the 
amplitude squared for the partonic process gp — > ggp, after averaging over the spin and 
color degrees of freedom, 

2 _27fa^Anfs 2 k 2 2 2 

lAl ~ — 16 — ^ ■ (84) 
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And the cross section for the partonic process gp — > ggp is 

da(gp^ggp) l f 2 2 27a 3 y k 2 2 2 1 

aivi x aK T 2 2 ^± 2 j |z s | , _ , (80j 



V ^ 

Following the same argument in the last subsection for the calculations of the partonic 
process qp — > qgp, we see that the dominant contribution of the integration over M\ comes 
from the region of M\ ~ 4k\^ where the differential gluon distribution function f g (x; Q 2 ) of 
the three terms in the integration X g can approximately take their values at the same scale 
of Q 2 = k\. That is, the integration X g is then simplified to 

1 1 3k 2 3 

T 9 = ^(-^(i + a k ))f g (x; k 2 T ) = Jjj^rf 9 ( x 'i k r) = Jpfg( x ^ k r)- ( 86 ) 



III. NUMERICAL RESULTS 

In this section, we give some numerical results for the diffractive light quark jet and 
gluon jet production at the Fermilab Tevatron. We will study the pt distribution and x± 
distribution of the cross section. A more thorough phenomenological study, including a 
comparison to currently available data at Tevatron on the diffractive dijet production rate, 
will be presented elsewhere. Again, we emphasize that our numerical results are just for the 
contribution from the "Pomeron Fragmentation" region (M| ~ 4k 2 -,). 

Another important effect which will affect the numerical results is the nonfactorization 
effect caused by the spectator interactions in the hard diffractive processes in hadron col- 
lisions mentioned in the section of Introduction. Here, we use a suppression factor Ts to 
describe this nonfactorization effect in the hard diffractive processes at hadron colliders |§. 
At the Tevatron, the value of Ts may be as small as Ts ~ 0.1 |5],|9]]- That is to say, the total 
cross section of the diffractive processes at the Tevatron may be reduced down by an order 
of magnitude due to this nonfactorization effect. In the following numerical calculations, we 
adopt this suppression factor value to evaluate the diffractive production rate of light quark 
jet at the Fermilab Tevatron. 

The cross section formulas for the three partonic processes have been calculated in the 
above section. We will use Eqs. (|24l) , (|58l) and fl85p to evaluate the production rates for 
these processes. In our calculations, the scales for the parton distribution functions and the 
running coupling constant are both set to be Q 2 = k\. For the parton distribution functions, 
we choose the GRV NLO set 131. 



The numerical results of the diffractive light quark jet production at the Fermilab Teva- 
tron are plotted in Fig. 5 and Fig. 6. In Fig. 5, we plot the differential cross section da/dt\t=o 
as a function of the lower bound of the transverse momentum of the light quark jet, A/rmm- 
This figure shows that the cross section is sensitive to the transverse momentum cut A/rmm- 
The differential cross section decreases over four orders of magnitude as /crmin increases from 
5 GeV to 15 GeV. 

It is interesting to compare the cross section of the diffractive light quark jet production 
with that of the diffractive charm quark jet production |13|], which is also shown in Fig. 5. 
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The two curves in this figure show that the production rates of the light quark jet and the 
charm quark jet in the diffractive processes are in the same order of magnitude. However, 
we know that the cross section of diffractive heavy quark jet production is related to the in- 
tegrated gluon distribution function, while the cross section of the light quark jet production 
is related to the differential gluon distribution function. 

In Fig. 6, we plot the differential cross section da/dt\ t =o as a function of the lower bound 
of the momentum fraction of the proton carried by the incident gluon £i m i n , where we set 
^Tmm = 5 GeV. This figure shows that the dominant contribution comes from the region of 
X\ ~ 10~ 2 — 10 _1 . This property is the same as that of the diffractive charm jet production 
at the Tevatron [JT3| . 

In Fig. 7, we plot the differential cross section da/dt\ t =o as a function of the lower bound 
of the transverse momentum of the gluon jet, £>rmin- This figure shows that the cross section 
is sensitive to the transverse momentum cut k Train . We plot separately the contributions 
from the two subprocesses, qp — > qgp and gp — > ggp. By comparison, we also plot the 
cross section of the diffractive light quark jet. The three curves in this figure show that 
the contribution from the subprocess gp — > ggp is two orders of magnitude larger than 
that from the subprocess qp — > qgp for the diffractive gluon jet production, and the light 
quark jet production rate is in the same order with that of the subprocess qp — > qgp. This 
indicates that the diffractive dijet production at hadron colliders dominantly comes from 
the subprocess gp — ► ggp in the two-gluon exchange model. 

In Fig. 8, we plot the differential cross section da/dt\t=o as a function of the lower bound 
of the momentum fraction of the proton carried by the incident gluon Xi m i n , where we set 
&Tmm = 5 GeV. Fig. 5 (a) is for the contribution from the subprocess qp — > qgp, and Fig. 5(b) 
is from the subprocess gp — > ggp. These two figures show that the dominant contribution 
comes from the region of X\ ~ 10~ 2 — 10 _1 for the subprocess gp — > ggp, and x\ > 1CT 1 for 
the subprocess qp — > qgp. These properties are similar to those of the diffractive charm jet 
and W boson productions calculated in [I3yi3|. 



IV. CONCLUSIONS 

In this paper, we have calculated the diffractive light quark jet and gluon jet productions 
at hadron colliders in perturbative QCD by using the two-gluon exchange model. We find 
that the production cross section is related to the squared differential gluon distribution 
function dG(x; Q 2 ) / dlnQ 2 at the scale of Q 2 ~ k\, where k? is the transverse momentum of 
the final state quark or gluon jet. We have also estimated the production rates of the light 
quark jet and gluon jet in the diffractive processes at the Fermilab Tevatron. 

As we know, the large transverse momentum dijet production in the diffractive processes 
at hadron colliders is important to study the diffractive mechanism and the nature of the 
Pomeron. The CDF collaboration at the Fermilab Tevatron have reported some results on 
this process [f|. Up to now, we have calculated all of the dijet production subprocesses in 
the diffractive processes at hadron colliders, including gp — * qqp, qp — > qgp and gp — ► ggp 
processes. In a forthcoming paper, we will compare the available data on the diffractive 
dijet production cross sections at the Tevatron |Sj with the predictions of our model to test 
the validity of perturbative QCD description of the diffractive processes at hadron colliders. 
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Figure Captions 



Fig.l. Sketch diagram for the diffractive dijet production at hadron colliders in perturbative 
QCD. The final state jet lines represent the outgoing quark or gluon lines, and the incident 
parton from the upper proton (labeled by Xipi) can be quark or gluon correspondingly. 
Fig. 2. The lowest order perturbative QCD diagrams for partonic process gp — > qqp. 
Fig. 3. The lowest order perturbative QCD diagrams for partonic process qp — > qgp. 
Fig. 4. The lowest order perturbative QCD diagrams for partonic process gp — > ggp. 
Fig. 5. The differential cross section da/dt\ t=0 for the light quark jet production in the 
diffractive processes as a function of kxmm at the Fermilab Tevatron, where ^ is the 
lower bound of the transverse momentum of the out going quark jet. For the charm quark 
jet production cross section formula, we take from Ref. and set m c = 1.5 GeV. 
Fig. 6. The differential cross section da/dt\ t= o for the light quark jet production as a function 
of £imin, where Xi min is the lower bound of Xi in the integration of the cross section. 
Fig. 7. The differential cross section da/dt\ t= o for the gluon jet production in the diffractive 
processes as a function of /cTmin- 

Fig. 8. The differential cross section da/dt\ t =o for the gluon jet production as a function of 
37-i m in, where xi m i n is the lower bound of x± in the integration of the cross section, (a) is for 
the contribution from the subprocess qp — > qgp, and (b) is from the subprocess gp — > ggp. 
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